Electrodiffusion of Lipids on Membrane Surfaces 
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Lateral translocation of lipids and proteins is a universal process on membrane surfaces. Local aggre- 
gation or organization of lipids and proteins can be induced when the random lateral motion is mediated 
by the electrostatic interactions and membrane curvature. Although the lateral diffusion rates of lipids 
on membranes of various compositions are measured and the electrostatic free energies of predetermined 
protein-membrane-lipid systems can be computed, the process of the aggregation and the evolution to the 
electrostatically favorable states remain largely undetermined. Here we propose an electrodiffusion model, 
based on the variational principle of the free energy functional, for the self-consistent lateral drift-diffusion 
of multiple species of charged lipids on membrane surfaces. Finite sizes of lipids are modeled to enforce 
the geometrical constraint of the lipid concentration on membrane surfaces. A surface finite element 
method is developed to appropriate the Laplace-Beltrami operators in the partial differential equations 
(PDEs) of the model. Our model properly describes the saturation of lipids on membrane surfaces, and 
correctly predicts that the MARCKS peptide can consistently sequester three multivalent phosphatidyli- 
nositol 4,5-bisphosphate (PIP2) lipids through its basic amino acid residues, regardless of a wide range of 
the percentage of monovalent phosphatidylserine (PS) in the membrane. 



INTRODUCTION 

Lateral motion of orientated lipids is an intrinsic mem- 
brane property. Through specific spatial and temporal 
variations of their lateral distribution the lipids regulate 
many essential membrane -related physiological processes. 
Multivalent phosphatidylinositol 4,5-bisphosphate (PIP2), 
for instance, plays a central role in anchoring proteins to 
plasma membranes and in regulating many ion channels, 
enzymes and cytoskeletons [1-4]. Charged lipids are also 
involved in driving membrane curvature that is required 
by cell endocytosis and exocytosis, as illustrated by the 
lateral movement of monovalent dioleoylphosphatidylser- 
ine (DOPS) to favor electrostatic interactions between neg- 
atively charged high-curvature membranes and positively 
charged concave surfaces of Bin/amphiphysin/Rvs (BAR) 
proteins [5, 6]. On the other hand, the lateral distribution 
of lipids are strongly affected by the interactions among 
lipid molecules, binding proteins, membrane curvature and 
membrane potential, and the ionic distributions in the solu- 
tion [5, 7-10]. It is therefore of great importance to study 
the lateral distribution of lipids under a wide range of phys- 
iological conditions to provide qualitative or quantitative 
information of related physiological processes. 

Lateral diffusion is a well recognized model of in-plane 
lipid motion [11-13]. Various techniques such as fluores- 
cence recovery after photobleaching (FRAP), nuclear mag- 
netic resonance (NMR), electron spin resonance (ESR), 
and fluorescence spectroscopy have been developed to de- 
termine the lateral diffusion coefficients of lipids [14-17], 
mostly by fitting the measured lengths of random walks to 
the linear dependence relation = Dt of normal 

diffusion to determine the the diffusion coefficient D. Sur- 
prisingly, the measured coefficients can differ by 2 orders 
of magnitude depending on the particular techniques, ma- 
terials and environments (in vitro or in vivo, for example) 



in experiments [18-20]. Furthermore, diffusion in mem- 
branes in vivo may not follow the normal linear relation be- 
cause the membrane contains a number of non-lipid com- 
ponents of different dimensions, which make the diffusion 
abnormal. According to the classical free-volume diffusion 
theory of Cohen and Tumbull [2 1-23], lateral diffusion in- 
deed follows a three-step procedure: (i) creation of a lo- 
cal free volume due to spontaneous thermal fluctuations; 
(ii) hopping of a diffusive molecule into the void; and (iii) 
filling of the hole left behind the hopping molecule by an- 
other molecule. Along with other attempts to describe the 
lipid diffusion in real complex membrane systems [_ i ], this 
theory has been integrated with the scaled-particle theory 
[25] to model the dependence of the diffusion coefficient 
on the sizes of lipids and non-lipid components [ ], and 
the theoretical results matched the experiments better than 
the original free-volume diffusion theory. Recent molec- 
ular dynamics (MD) simulations of membrane dynamics 
have provided different insights of the mechanism of lat- 
eral diffusion of lipids. It was shown that the motion of 
neighboring lipids are strongly correlated, therefore mo- 
tion of lipids in clusters is more favored over the hopping 
of individual lipids [27]. 

Applications of these measured lateral diffusion coeffi- 
cients to computational simulations or quantitative analysis 
of physiological processes, however, are scarce, probably 
due to a lack of appropriate models of the mechanisms that 
can lead to the variation of the lateral distribution of lipids. 
Suppose that the lateral diffusion of lipids in a bounded ple- 
nary 2-D domain 17 is described by the standard diffusion 
equation with source 

^ = DAp + f{x,y), (1) 

where D is the constant diffusion coefficient, then the vari- 
ation of the continuous lipid concentration p can be caused 
by its initial non-uniform distribution, the product function 
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f{x, y) that models the addition or removal of lipids in the 
domain, or the prescribed distribution of p on the boundary 
of Vt. None of these three features have been well mod- 
eled for lipid diffusion. We note, however, that similar 
models are developed for lateral diffusion of membrane 
proteins. For instance the aggregation of proteins to the 
growing clusters in membranes is modeled as an absorbing 
boundary condition on the moving boundary [2;., |. For 
lateral diffusion of lipids or proteins on curved biological 
membranes one needs to replace the Laplace operator in 
Eq.(l) by the Laplace-Beltrami operator so the effects 
of surface curvatures on the diffusion rate can be incor- 
porated [30-32]. Such a model also enables us to inves- 
tigate the long-term effects of the random fluctuations of 
the surface metric on the lateral diffusion [33]. More im- 
portantly, membrane proteins and many regulatory lipids 
such as PIP2 and DOPS are charged hence a drift-diffusion 
equation is more suitable than the standard diffusion equa- 
tion (1). Electrodiffusion of charge particles in solution 
has been extensively studied in past decades, and there has 
been resurgent interests recently in modifying the classi- 
cal Poisson-Nernst-Planck (PNP) theory to better model 
the drift-diffusion of charged particles with finite sizes and 
particle -particle correlations [34, 35]. PNP theory and its 
modifications have been extended recently to study the lat- 
eral motion of lipids or proteins in membranes. For ex- 
ample, 2-D Poisson-Boltzmann-Nernst-Planck equations 
with size exclusion are derived from a constraint free en- 
ergy minimization to study the lateral motion of lipids in 
the membrane [36], assuming the co-existence of mobile 
ions, peptides, and lipids in the hydrated leaflets of the 
membrane. Formation and disassociation of peptide-lipid 
complexes are modeled via reaction terms in the Nernst- 
Planck equations, and results are compared favorably to 
their dynamic Monte Carlo simulations [37]. 2-D PNP 
equations without size exclusion were used to simulate the 
electrodiffusion of lipids on a planar surface [38], coupled 
with a nonlinear 3-D Poisson-Boltzmann equation with 
a charge density neglecting the contribution of charged 



where pi is the ion concentration of the i*'' species and is 
the lipid concentration of j*'' species, with i = denotes 
the neutral solvent (water) molecules and j = denotes 
the neutral species of lipids. The characteristic function A 



lipids. While various degrees of agreement to the ex- 
perimental measurements or other theoretical predictions 
have been achieved with these models, an electrodiffu- 
sion model similar to the modified Poisson-Nemst-Planck 
equations in bulk are necessary for quantitative study of 
self-consistent distributions of charged particles on general 
curved surfaces. 

In this paper we developed a generic surface drift- 
diffusion equation based on a generalized Borukhov model 
[35, 39] to model the electrostatic mediated lateral motion 
of charged lipids on arbitrarily curved membrane surfaces. 
The bilayer membrane is described as a dielectric contin- 
uum with a constant permittivity and continuous distribu- 
tions of surface charges, the latter model the distributions 
of charged lipids on membrane surfaces. The lipids are 
treated as diffusive hard disks, and their effective radius 
are taken into account together with the finite sizes of mo- 
bile ions in the entropic contribution of the total free en- 
ergy. Without proper modeling of the lipid sizes the lat- 
eral concentration of charged lipids can be easily overesti- 
mated wherever there are strong attractive electrostatic in- 
teractions, even though the molecular dimensions of lipids 
have been considered in the determination of the lateral dif- 
fusion coefficient. We derive a system of nonlinear partial 
differential equations (PDEs) and develop the related nu- 
merical algorithms. Numerical simulations are carried out 
to study the surface electrodiffusion of PIP2 due to the elec- 
trostatic interactions between the membrane and the MAR- 
CKS protein under a wide range of physiological condi- 
tions. Our results are well validated by experimental mea- 
surements and atomistic modeling of lipid sequestration. 

MATHEMATICAL MODELS OF SURFACE 
ELECTRODIFFUSION 

We consider a solvated protein-membrane system as il- 
lustrated in Figure 1(A). In the spirit of the classical fluid 
density functional theory [35, 36, 40] we define the fol- 
lowing free energy for the membrane-protein-ionic solu- 
tion system: 



(2) 



I 

is 1 in the solution domain Vtg and zero elsewhere. X is the 
coordinates of the membrane surfaces St and S]^. kB,T,(j) 
are the Boltzmann constant, temperature, and the electro- 
static potential, respectively. Zi , are the valences of cor- 
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responding ions and lipids, and e is the elementary charge, 
is the permanent charges in proteins. 

1 

is the fugacity of the i^^ species of ions with effective size 
Oi, and 



1 



is the fugacity of j*'* species of lipids approximated as hard 



Oo and aj,, and their concentrations poi Po follow the rela- 
tions 



pQ<A + ^P^o!l = 1, 

1=1 



(3) 
(4) 



These two constraints are enforced directly in the varia- 
tions of F with respect to pi and p[ without resorting to a 
Lagrangian multiplier, in contrast to [ ^ f ], giving rise to the 



disks with an effective size [35, 41]. Here the solvent electrochemical potentials p^ for ions (i > 1) and p - for 
molecules and neutral lipids have respective effective sizes charged lipids (j > 1): 
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where = (flj/ao)^, A;'- = [a^Ja^Y. These electro- 
chemical potentials are related to the flux of ions and lipids 
through the constitutive relations for fluxes 
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where is the surface gradient operator, and irij are 
the mobilities of respective species of ions and lipids; they 
are related to the diffusivities Di or Dj through Einstein's 
relation 

Di = rriikBT or = m[jkBT. 

The drift-diffusion equations for ions are obtained from the 
mass conservation in solution: 
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and the surface drift-diffusion equations for lipids are ob- 
tained from the mass conservation on the membrane sur- 
faces if it evolves at a given velocity field Ug [42]: 
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where n is the outer normal vector of the surface. If there 
are only two types of lipids of the same effective size, one 
charged and the other neutral, in a macroscopically static 
membrane, one can solve only one surface drift-diffusion 
equation, which is now simplified to be 

(9) 

in which the subscription j for species is no longer needed. 
For enclosed membrane surfaces Equations (8) and (9) do 
not have boundary conditions, but both are subject to the 
constraint of mass conservation 

/ Pjds = T,, (10) 
Js 

where T, is the given total quantity of a species of charged 
lipids on the surface S. Extremization of F with respect to 
(f) using variation 6F/6(p gives the Poisson equation for the 



electrostatic potential (j): 

-\/-{eV(j)) = pf + X^p,z,e inn, (11) 

i=l 

with interface conditions 

ernf,-^^ = ^s-^ + 2^ Pj^je> 4>s = (Prnb OU StUSb, 
J = l 

e™^— = ^s^— 9s = 4>m on i , 
on on 

where the charge density carried by charged lipids on the 
membrane surfaces is written as an interface condition. 
The subscript m, m6, s denote the quantities in proteins, 
membrane, and the solution, respectively. 

To acquire an accurate characterization of the diffusion 
process we will solve the Poisson equation (11) coupled 
with the time-dependent drift-diffusion equation (7) and 
the time-dependent surface drift-diffusion equation (8), us- 
ing a finite element method [3 ] and a linear surface finite 
element method [43]. If only the equilibrium distributions 
of ions and charged lipids are of interest one can instead 
solve the Poisson equation coupled with steady-state drift- 
diffusion equation and stead-state surface electrodiffusion 
equation. Notice the latter. 
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admits a trivial solution = if the constraint (10) is not 
enforced. To obtain a physical solution that is consistent 
with this constraint we introduce the decomposition 



where 
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is the known average concentration of p'- on the surface S 
with total surface area IS*!. Replacing p'- in equation (12) 
with this decomposition we get 
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from which the nontrivial variation and hence the non- SOLVATED MARCKS-MEMBRANE SYSTEM 

trivial p^ can be uniquely solved. 

The system, as illustrated by Figure l(A-B), consists 
of a MARCKS peptide and a membrane vesicle whose 
interior and exterior radii are 230A and 270A, respec- 
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FIG. 1. Solvated MARCKS-membrane system and the com- 
putational model. (A) Schematic illustration of the MARCKS 
peptide residing on the inner leaflet of a vesicle. Lipids with 
charged (red) and polar (cyan) headgroups are distributed in the 
membrane. (B) 2D cross section of the computational model 
for (A). The domain f2„jf, is the bilayer membrane modeled as 
a dielectric continuum without atomistic details. Distribution of 
lipids on membrane surfaces St and Si, follows the surface elec- 
trodiffusion equation. Mobile ions are distributed in the aqueous 
solution exterior and interior to the vesicle. Atomistic structure 
of the MARCKS peptide is retained and singular charges are dis- 
tributed in Qm- r is the molecuar surface of the peptide, and 
dQ is the exterior boundary of the computational domain for the 
Poisson equation. (C) 3D tetrahedral discretization of the entire 
computational domain. (D) local refinement of the tetrahedral 
mesh near the MARCKS peptide. 

lively. The MARCKS peptide is unfolded and the pep- 
tide is built from the sequence of the effective domain 
as given in [44]. The peptide is placed above the in- 
ner leaflet of the vesicle, the closest distance between the 
peptide and the membrane surface is about 2A, with the 
assumption that the position and orientation are not af- 
fected by its interaction with the membrane and lipids. 
The insertion of MARCKS is not modeled here because 
the aromatic phenylalanine residues penetrating to the acyl 
chain region obstruct the diffusion of lipids. Experi- 
ments and Monte Carlo simulations suggest that the dif- 
fusion coefficients measured on systems free of inserted 
molecules may not be applicable to the obstructed diffu- 
sion [45, 46]. In order to investigate electrodiffusion of 
lipids with different valences we will consider two types 
of lipid compositions, PIP2/PC(Phosphatidylcholine) and 
PIP2/PS(phosphatidylserine)/PC; all head groups are as- 
sumed to have the same effective diameter 8.3A. The two 
leaflets of the membrane are assigned the same average 



lipid composition. The dielectric constant = e^b = 2 
in the peptide and the membrane, and €s = 80 in the solu- 
tion. There are significant disparities in the measured lat- 
eral diffusion coefficients of lipids in the plasma membrane 
[15, 47], here we choose = 580 / jis. The solution 
contains the salt KCl only, whose concentration will be ad- 
justed between 50mM to 300mM for the examination of the 
effects of the salt concentration on the lateral electrodiffu- 
sion. The effective sizes of K and CI are 2.7Aand3.6A, re- 
spectively, and their diffusion coefficients are 78000A^/;U. 
We note that the final equilibrium distributions of lipids 
and mobile ions are independent of these diffusion coeffi- 
cients. There is no ion exchange through the membrane. 
The Poisson equation (11) is solved in the entire 3-D do- 
main Vis U r^nif, U rj,„, the bulk drift-diffusion equation (7) 
is solved in the 3-D solvent domain fi^, and the surface 
drift-diffusion equation is solved on both spherical mem- 
brane surfaces St and 5f,. This way the lateral diffusion of 
lipids on two membrane leaflets are coupled. A cross sec- 
tion of the finite element tetrahedral mesh and the locally 
refined mesh near the peptide are shown in Figure l(C-D). 



COMPUTATIONAL SIMULATIONS AND DISCUSSIONS 

Effects of finite size of PIP2 on its surface distribution 

We compare the aggregation of two models of PIP2 
lipids due to the electrostatic attraction of the MARCKS 
peptide. In one model the lipids are described as parti- 
cles with vanishing sizes and in the other the lipids have 
an effective diameter 8.33A. We observe the aggregation 
by starting with an initial uniform distribution of 1 % PIP2 
on membrane surfaces. The histories of aggregation are 
shown in Figure 2. The accumulation of charged lipids 
is found much faster when the size effects are neglected, 
and its concentration will quickly exceed the upper limit 
of the lipid concentration, resulting in an unphysical dis- 
tribution of lipids on membrane surfaces. This upper limit 
of the concentration of PIP2 (i.e., the membrane consisting 
of 100% PIP2) is 0.0144/A2, corresponding to about 36 
lipids on a 50x50A^ membrane surface [5]. The overesti- 
mated concentrations of charged lipids can lead to an over- 
estimated dielectric surface force density in studying the 
protein-membrane electrostatic interactions [48]. In con- 
trast the model with lipid size correction predicts a lipid 
concentration that is well bounded by the saturation value. 



Dependence of PIP2 aggregation on its average 
concentration 

Experiments [1, 2, 4] show that a MARCKS peptide 
can consistently sequester three or four PIP2 lipids when 
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FIG. 2. Without an appropriate model of the lipid size the 
PIP2 quickly exceeds the upper limit of the concentration 
under the attraction of the MARCKS peptide. Blue: lipids 
with vanishing size; Green: lipids with an effective diameter 
8.33A. The red line signifies the limit of the lipid concentration 
0.0144/A^. Initial average percentage of PIP2 is 1%. The scale 
on the vertical axis is amplified 100-fold. 

the average fraction of PIP2 in the membrane varies be- 
tween 0.01% and 1%. Suppose 1% PIP2 is uniformly dis- 
tributed in the membrane, 3 PIP2 lipids shall occupy a sur- 
face area of about 144 x 144A^. Such an area is nearly 33 
times larger than the area coverage of MARCKS peptide 
(~45 X ISA^). For 0.1% PIP2 this area amounts to almost 
1/3 of the total area of the inner leaflet of the vesicle with 
an inner radius 230A; and there is barely one PIP2 lipid 
on the entire inner surface of the vesicle when the mem- 
brane contains only 0.01% PIP2- Consequently we define 
the lower bound of the PIP2 percentage in our simulations 
to be 0.1%. The maximum percentage of PIP2 in our sim- 
ulations is 30%, for which 3 PIP2 lipids occupy an area 
of 26.4x26.4A^ when distributed uniformly; this area ap- 
proximates the coverage of MARCKS peptide on the mem- 
brane surface. This general assessment of the sparsity of 
PIP2 lipids suggests that the equilibration of PIP2 distribu- 
tion in response to the electrostatic perturbation of MAR- 
CKS peptide will vary as the initial percentage of PIP2 
changes. To quantify these differences and to reveal the 
continuous dependence of the sequestration of PIP2 on its 
average concentration we simulate the electrodiffusion of 
PIP2 with five initial conditions; the results are summa- 
rized in Figure 3. With 30% PIP2 it takes about 0.2ps to 
re-equilibrate its distribution in response to the approach- 
ing of MARCKS peptide, while it takes about 10/is if the 
membrane contains 1 % PIP2 . For an initial percentage of 
0.1% this time becomes as long as 100 fis. At a low per- 
centage of PIP2 it takes a much longer time for the lipids 
to move over a larger distance to the vicinity of the MAR- 
CKS peptide so the lipids can be sequestered. Cautions 
must be taken in interpreting this long time of aggregat- 
ing because the attraction of PIP2 to the peptide is not 
always started with a uniform distribution of PIP2. In- 
deed, the MARCKS peptides are recycled in real biolog- 
ical systems: when [Ca^+] increase in the local cytoplasm 
the Ca2+/calmodulin (Ca/CaM) will bind to the MARCKS 
peptides that have attached to the membrane and pull them 



off. When the detached MARCKS proteins move away 
from the membrane the sequestered PIP2 lipids will be re- 
leased and diffuse away from the aggregated region [8]. 
If the period of this recycle is sufficiently small, the re- 
attached MARCKS peptide will quickly attract the PIP2 
lipids back before the re-establishment of an uniform distri- 
bution of PIP2. It is therefore possible to integrate the cur- 
rent continuum surface electrodiffusion equations for lipids 
with the diffusion-reaction dynamical models for Ca^^, 
calmodulins, and other associated molecules to yield an ac- 
curate characterization of the biochemical signals that are 
transmitted through these processes. 

We count the number of sequestered PIP2 lipids by inte- 
grating the density on the membrane surface 

^PiP2 = / P^ds, 

JSp 

where the sequestration domain Sp is defined to be the 
projection to the inner membrane surface of a 45xl5A^ 
rectangular centered at Serl62:N, as shown in Figure 3(C). 
This rectangular is chosen to enclose the entire MARCKS 
peptide. We note that this definition of sequestration is 
different from those in the atomistic modeling of lipids 
[1, 49], where PIP2 lipids can be represented as potential 
surfaces of -25mV. The sequestration of lipids is identi- 
fied as the enclosure of these negative potential surfaces 
by the iso-surfaces of positive potential (~25mV) that are 
induced by the basic residues in the peptide. Quantitative 
agreements are found between the current continuum and 
the previous atomistic models of lipids. At 0.1% initial 
percentage, about 2 lipids will finally be attracted to the 
membrane surface under the peptide after a long travel, 
and about 3 PIP2 lipids can be sequestered at 1% initial 
percentage. More PIP2 lipids are found to be clustered be- 
low the peptide with a further increase of its percentage, 
gradually saturating to the maximum density of lipids. Fig- 
ure 3(B) indicates that the maximum concentration of PIP2 
reaches the saturation density at 30% of average concen- 
tration. Correspondingly there are about 7 PIP2 are clus- 
tered below the peptide, c.f. Figure 3(D), which account for 
72% of the total lipids in the domain Sp. It is worth noting 
that this aggregated distribution of PIP2 is generated from 
the minimization of the total free energy (2), and is consis- 
tent with the electrostatic potential field of the membrane- 
MARCKS system by construction. In contrast, with atom- 
istic models of the sequestration one needs to place charged 
lipids at positions that are determined a priori so that the 
electrostatic potential can show favorable electrostatic in- 
teractions between basic peptides and phospholipid mem- 
branes [:: ]. The current continuum model thus has a great 
potential in modeling the membrane -protein systems for 
which the positions of specific lipids are unknown a pri- 
ori or change dynamically with the conformational change 
of membrane or proteins during their interactions. 
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FIG. 3. (A) Accumulation of PIP2 near the MARCKS pep- 
tide; (B) the dependence of the maximum concentration of 
PIP2 on the initial average concentration; (C) Integration re- 
gion used for counting the sequestered PIP2; and (D) Number 
of sequestered PIP2 . The initial average concentrations of PIP2 
in (A) are 1%, 5%, 10%, 20%, 30% for the green, blue, black, 
cyan, and maroon lines, respectively. Complete history for the 
simulation with initial 1% PIP2 is shown in Figure 2. In (C) the 
red lines represents the MARCKS peptide above the membrane 
surface. The integration region is not pre-determined for a con- 
forming triangulation so its boundary is not straight. The average 
percentages in (D) are 0.1%, 1%, 5%, 10%, 20% and 30%. 

Dependence of PIP2 aggregation on the concentration of 
monovalent lipids 

When there are multiple species of negatively charged 
lipids on the membrane their concentrations near the MAR- 
CKS peptide can be completely different. Boltzmann re- 
lation predicts that multivalent lipids such as PIP2 have 



an affinity that is 1000 times stronger than the monova- 
lent lipids such as phosphatidylserine (PS) when binding 
to the MARCKS peptide. Although the real distribution of 
charged lipids on a membrane does not follow the Boltz- 
mann relation due to their large size and correlations, elec- 
trostatic free energy of the sequestration showed that PIP2 
is favored by the peptide compared to PS [31]. To quan- 
tify this competition we simulate electrodiffusion of PIP2 
and PS with different fractions. The average percentage of 
PIP2 is fixed at 1 % and the fraction of PS varies between 
1% and 30%. Figure 4 shows the weak dependence of the 
aggregation of PIP2 on the fractions of PS that are under in- 
vestigation. While the total concentrations of PIP2 and PS 
increase to approach the saturation value with the increase 
of average concentration of PS, the peak concentration of 
PIP2 is found decreasing. This decrease is slight never- 
theless. The number of sequestered PIP2, computed using 
equation (1), reaches it smallest value of 2.7 when the aver- 
age percentage of PS is 30%; the corresponding number of 
sequestered PS lipids is less than 0.4. These results are con- 
sistent with the experimental observations that monovalent 
acidic lipids are less likely be sequestered by membrane- 
bound basic peptides if there are multivalent acidic lipids 
on the membrane [51], and that the sequestration of PIP2 
is more favorable when the mol percentage of monovalent 
acidic lipids in the membrane decreases [50]. 

The time dependent solutions of electrodiffusion pro- 
vides quantitative information that can not be attained 
through the investigation of electrostatic free energy only. 
An examination of Figures 4(C-E) shows that there is an 
increase of monovalent PS lipids for the first few microsec- 
onds, in contrast to the continuous growth of the concentra- 
tion of multivalent PIP2. The strength of positive electro- 
static potential near the sequestration region is reduced due 
to the accumulation of negatively charged lipids. It appears 
that the monovalent PS lipids are more susceptible to this 
local change potential so their concentration starts to de- 
crease after a maximum value is reached. This temporary 
increase of monovalent PS lipids is negligibly small when 
their average percentage is low, c.f. Figure 4(A). 

Dependence of PIP2 aggregation on the ion concentration 

Experiment observations suggest that the aggregation 
of charged lipids depends on the local electrostatic field 
[49, 52], whose strength can change drastically with the 
ion concentration in the solution. This dependence is mod- 
eled in the current work through the coupling of the surface 
electrodiffusion equation with the Poisson equation (11) 
that admits the ion concentrations determined by the mod- 
ified Nernst-PIanck equations (7). To validate the model 
we simulate the PIP2 aggregation with KCI concentration 
ranging from 50mM to 300mM. Figure 5 illustrates that 
PIP2 lipids are less aggregated for large salt concentrations, 
as a result of the electrostatic screening and the weakened 
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FIG. 4. (A-E) Accumulation of PIP2 and PS near the MAR- 
CKS peptide and (F) the dependence of the maximum con- 
centrations of PIP2 and PS on the initial average concentra- 
tion of PS. The initial average concentrations of PS in (A-E) are 
1%, 5%, 10%, 20%, 30% respectively, with blue lines denoting 
PIP2 and red lines denoting PS. The scale on the vertical axis is 
amplified 100-fold. 
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FIG. 5. Salt-dependent lipid sequestration. Maximum 
concentration of PIP2 decreases as salt concentration increases 
from 50mM (green), lOOmM (magenta), I40mM (blue), 200mM 
(cyan) to 300mM (red). The number of sequestered PIP2 lipids 
are 4.02, 3.35, 2.89, 2.16 and 1.07, respectively. 



with the experimental observations and atomistic model- 
ing of lipid sequestration. The ability to efficiently and ac- 
curately identify the local concentrations of specified pro- 
teins or lipid headgroups will enable us to integrate the self- 
consistent lateral translocation of lipids and proteins to the 
analysis of membrane trafficking and curvature, and related 
signal pathways and diseases. 



tangent component V^f/i) of the electrostatic field. Our sim- 
ulations also show that other species of mobile ions such 
as Ca^^, Mg^^ will have similar screening effects. Note 
that the Ca/CaM can also reverse the binding of MAR- 
CKS peptide to the membrane, so the current model needs 
to be integrated with models for Ca/CaM association and 
Ca/CaM/MARCKS binding [53] to give a full description 
of the effects of Ca^+ on PIP2 sequestration. 



SUMMARY 

Lipid headgroups and membrane proteins diffuse later- 
ally on the membrane surfaces. Local aggregation of lipids 
and proteins is usually a precursor of membrane fusion, 
fission, and budding and therefore is a critical step in many 
signaling pathways that are regulated by the membrane cur- 
vature. Lateral diffusion of charged headgroups are medi- 
ated by the nonspecific electrostatic interactions and hence 
can not be model by the conventional diffusion equation. 
Here we establish a surface electrodiffusion model with 
finite-size correction to describe the self-consistent lateral 
aggregation of charged lipid headgroups on curved mem- 
brane surfaces due to membrane/protein electrostatic inter- 
actions, without requiring the knowledge of locations of 
individual aggregated lipids. Our results agree very well 
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